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1. Introduction 



The complex quaternion algebra (biquaternion algebra) Q is well known as a four dimensional 
vector space over the complex number field C with its basis 1, ei, 62, 63 satisfying the muhiph- 
cation laws 



= 62 = 63 = -1, 616263 = -1. 

6162 = -6261 = 63, 6263 = -6362 = 61, 6361 = -6363 = 62, 

and 1 acting as unity element. In that case, any element in Q can be written as 

a = ao + aiei + 0262 + 0363, 



(1.1) 
(1.2) 

(1.3) 



where bq — 03 G C. According to this definition, real numbers, complex numbers, and real 
quaternions all can be regarded as the special cases of complex quaternions. A well-known fun- 
damental fact on the complex quaternion algebra Q (see, e. g., [@, |6|, 0]) is that it is algebraically 
isomorphic to the 2x2 total matrix algebra C^^^ through the bijective map ^p : Q — > C^^^ 
satisfying 



1 
1 



V'(ei) 







V'(e2) 



-1 







'0(63) 











These four matrices are well-known as Pauli matrices. Based on this map, every element a 
ao + aid + 0262 + 0363 S Q has a faithful complex matrix representation as follows 



ao + ail 
a2 - a^i 



-(02 + asi] 
ao — aii 



G 



^2x2 



(1.4) 
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In this article, we shall reveal a deeper relationship between a and V'(o), which can simply be 
stated that there is an independent invertible matrix Q of size 2 over Q such that all a G Q 
satisfy the following universal similarity factorization equality 

Q^^(i\a.g{a, a)Q = ipia), 

where Q has no relation with the expression of a. Moreover we also extend this equality to all 
mxn matrices over Q. On the basis of these results, we shall consider several basic problems re- 
lated to complex quaternion matrices, such as, generalized inverse, eigenvalues and eigenvectors, 
similarity, and determinant of complex quaternion matrices. 

Some known terminology on complex quaternions are listed below (see, e.g., Q). For a = 
ao + aiei + 0262 + 0363 € Q, the dual quaternion of a is 

a = ao - fliei - 0262 - 0363; (1.5) 

the complex conjugate of a is 

a* = + olei + 0^62 + 0363; (1-6) 

the Hermitian conjugate of a is 

= (a)* = a^J- olei - 0^62 - 0^63; (1.7) 

the weak norm of a is 

n{a) = aQ + af + a2 + a^. (1.8) 

A quaternion a S Q is said to be real if a* = a, to be pure imaginary if a* = —a, to be scalar if 
a = a, to be Hermitian if = a. 

For any A = (a^j) G Q™^", the dual of ^ is ^ = (otj) G Q"^™-; the Hermitian conjugate of 
A is = (alg) S Q"^™". A square matrix A is said to be self-dual if ^ = ^, it is Hermitian if 
j4"f = A, it is unitary if AA^ = A'^A = I, the indentity matrix, it is invertible if there is a matrix 
B over Q such that AB = BA = I. 

Some known basic properties on complex quaternions and matrices of complex quaternions 
are listed below. 

Lemma l.l|l|]|^]. Let a, b £ Q be given. Then 
(a) a = a, (a*)* = a, (a^)^ = a. 

(h) ^Tb = a + b, (a + b)* = a* + b*, (a + 6)t = ot + 6t. 

(c) ^ = ba, (ab)* = a*b*, (o6)t = b^al 

(d) aa = aa = n{a) = n{a); 

(e) a is invertible if and only if n{a) 7^ 0, in that case a^^ = n^^{a)a. 

Lemma 1.2§. Let A G Q"^^", B G Q"^p be given. Then 

(a) J = A, (At)t = A. 

(b) AB = KA, {AB)^ = B^Al 
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(c) (AB) ^ = B ^, if A and B are invertible. 

(d) (A)-^ = (]FT), (^t)-i = (A-i)t, if A is invertible. 

2. A universal similarity factorization equality over complex 
quaternion algebra 

We first present a general result on the universal similarity factorization of elements over 
2x2 total matrix algebra. 

Lemma 2.1. Let M2(F) be the 2x2 total matrix algebra over an arbitrary field F with its 
basis eii, eu, 621 and 622 satisfying the following multiplication rules 



-sq, 



p 



0, t^p 



s, t,p,q = 1, 2. 



(2.1) 



Then for any a = ancn + 012612 + 021621 + 022622 & M2{T), where Ost G ^, the corresponding 
diagonal matrix diag( a, a) satisfies the following universal similarity factorization equality 



Q 



a 
a 



where Q has the independent form 



Q = Q-^ = 



Oil ai2 
021 022 



611 621 

ei2 622 



•2x2 



(2.2) 



(2.3) 



Proof. According to Eq.(2.1), it is easy to verify that the unity element in M2{J-) is 
6 = 611 + 622- 111 that case, the matrix Q in Eq.(2.3) satisfies 



611 


621 




. 612 


622 . 





611 + 621612 611621 + 621622 
. 612611 + 622612 612621 + 622 



611 + 622 
611 + 622 



e/2, 



which implies that Q is invertible over M2 (J-) and Q = Q ^ . Next multiplying the three matrices 
in the left-hand side of Eq.(2.2) yields the right-hand side of Eq.(2.2). □ 

Theorem 2.2. Let a = oq + ai6i + 0262 + 0363 £ Q be given. Then the diagonal matrix 
diag( a, a) satisfies the following universal factorization similarity equality 



Q 



a 
a 



Q 



aQ + aii -(02-1-03^) 



2x2 



(2.4) 



where Q is an unitary matrix over ' 



1 - i6i 62 + i63 

—62 + ^63 1 -I- zei 



(2.5) 



Proof. According to Eq.(2.4), we choose a new basis for Q as follows 

611 = ( 1 - ■iei ) , ei2 = ;^ ( -62 + ^e3 ) , 621 = ( 62 zes ) , 622 = ( 1 + ^62 ) . (2.6) 
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Then it is not difficult to verify that the above basis satisfies the multipUcation rules in Eq.(2.1). 
Under this new basis, any element a = oq + oiei + 0262 + 0363 G Q can be expressed as 



a = {ao + iai )eii + ( -02 — ia^ )ei2 + ( 02 - ias )e2i + ( ao - iai )e22- 
Substituting Eqs.(2.6) and (2.7) into Eq.(2.2), we obtain Eqs.(2.4) and (2.5). □ 
The equality in Eq.(2.4) can also equivalently be expressed as 



(2.7) 



Q 



an 


ai2 


Q = 


a 


" 






_ 




. 0,21 


(122 . 


a 



v>2x2 



where Ggt G C is arbitrary, Q is as in Eq.(2.5), and a has the form 

a = ^(aii + «22 ) + ^("22 - an )iei + ^(«2i - 012 )e2 + ^(^12 + 021 )ie3. 



(2.8) 



(2.9) 



This equality shows that every 2x2 complex matrix is uniformly similar to an diagonal matrix 
with the form a/2 over the complex quaternion algebra Q. 

The complex quaternions and their complex matrix representations satisfy the following 
operation properties. 

Theorem 2.3 Let a = qq + aici + 0262 + a3e3, 6 G Q, A G C be given. Then 

(a) a = b <J=^ V(a) = V'(&)- 

(b) ^(a + 6) = V;(a)+V(6), ^(ab) = ^{a)^{b), V(Aa) = V(aA) = AV'(a), = h- 



(c) V'(a) = 

(d) ^(a*) = 



1 
-1 

1 

-1 



r («) 



-1 

1 

-1 

1 



(e) ip{a^) = i^{a) = ip*ia), the conjugate transpose oftp{a). 

(f) det ipla) = n{a) = ag + af + 03 + ag; 

(g) a = \E2'4}{a) e\, where E2 = [I - iei, 62 + ie^]. 

(h) a is invertible if and only if 'ip{a) is invertible, in that case, ip{a~^) = ip~^(a) and 
a-^ = lE2ij-Ha)El 

For a noninvertible element over Q, we can define its Moore-Penrose inverse as follows. 
Definition. Let a G Q be given. If the following four equations 

(ax)^ = ax, (xa)^ = xa (2-10) 



axa = a, 



have a common solution x, then this solution is called the Moore-Penrose inverse of a, and 

denoted by x = a+. 

The existence and the uniqueness of the Moore-Penrose inverse of a complex quaternion a 
can be determined by its matrix representation ^(a) over C. In fact, according to Theorem 
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2.3(a), (b) and (e), the four equations in Eq.(2.10) are equivalent to the fohowing four equations 
over C 

ip{a)ip{x)'ip{a) =i/j{a), tlj{x)ip{a)'ip{x) =i/j{x), 
[^(a)^(x)]* = ^(a)^(x), [^(x)^(a)]* = ^(x)^(a). 
According to the complex matrix theory, the following four equations 

^{a)Y^P{a) = i/jia), Yil;{a)Y = Y, [Ha)Y]* = ^{a)Y, [Yi;{a)]* = Yil;{a). 

has a unique solution Y = V'^(a)) the Moore-Penrose inverse of tpi^')- Then by Eq.(2.8), it 
follows that there must be a unique x over Q such that tpi^) = Y = ilj'^{a), in which case, this 
X can be expressed as 

x = ^E2Y4 = ^E22p+ia)El. 

Correspondingly this x is the unique solution to Eq.(2.10). In summary, we have the following. 

Theorem 2.4. Let a E Q be given. Then its Moore-Penrose inverse a"*" exists uniquely, and 
satisfies the following four equalities 

V^(a+) = V+(a), a+ = V+W 4, 

where E2 = [1 — iei, 62 + ie^]. 

One of the basic problems related to complex quaternions is concerned with similarity of two 
complex quaternions. As usual, two complex quaternions a and b are said to be similar if there is 
an invertible complex quaternion x such that x~^ax = b, and this is written as a ~ 6. It is easy 
to verify that the similarity mentioned here is an equivalence relation on complex quaternions. 

A simple result follows immediately from the above definition, Eqs.(2.4) and (2.8). 

Theorem 2.5. Let a, b e Q be given. Then 

a ~ 6 <^ il){a) ~ ip{b). (2.11) 



Prom Eq.(2.11), we easily find the following. 

Theorem 2.6. Let a = + aici + 0262 + 0363 € Q given with a ^ C 

(a) // + + 03 / 0, then a ~ ao + r(a)ei, where T(a) is complex number satisfying 
T^(a) = af + 02 + «!• 

(b) // af + 02 + 03 = 0, then a 



~ ao - 562 + 5^63. 



Proof. Por any a G 

ijj{a) is 

I XI2 - ^{a) I = 



J, the characteristic polynomial of its complex matrix representation 

A - (oo + aii) 



-0.2 + a^i 



a2 + 031 
A — ( ao — aii ) 



(A-ao) +0^ + 02 + 03. 



5 



From it we immediately know that if af + 02 + a| 7^ 0, then 



ao + r(a)i 

ao — T{a)i 



iplao + r(a)ei ], 



(2.12) 



and if af + 03 + ag = 0, then 



ao 1 
ao 



1 1 • 

ip { ao - + 2^63 



(2.13) 



Correspondingly applying Eq.(2.11) to Eqs.(2.12) and (2.13) may lead to Part (a) and Part (b) 
of this theorem. □ 

3. Two universal factorization equalities on complex quaternion 
matrices 

In this section, wc extend the universal similarity factorization equality in (2.4) to any mxn 
matrix over Q, and give some of its consequences. 



Theorem 3.1. Let A = Ao + Aici + ^262 + ^363 G Q"'''" 
following universal factorization equality 



Q 



2m 



A 

A 



Q 



2n 



Ao + Aii -{A2 + A3i' 
A2 - Asi Ao - Aii 



be given. Then A satisfies the 
:= -^{A) e C2"»x2", (3.1) 



where Q2t has the independent form 

1 



Q2t = Q2t = ^2* - 2 



( 1 - iei )It ( 62 + ies )It 
( -62 + ies )It ( 1 + iei )It 



t = m, n. (3-2) 

In particular, when m = n, Eq.(3.1)becomes a universal similarity factorization equality over Q 



Proof. It follows directly from multiplying out the three block matrices in the left-hand side 
ofEq.(3.1). □ 

The matrix ^{A) in Eq.(3.1) is called the complex representation of A. If setting 
An = Ao + Aii, Ai2 = -{A2 + Asi), ^21 = ^2 - ^3^, An = Ao - Aii 
in Eq.(3.1), then it can equivalently be expressed as 



Q 



2m 



' Au 


A12 ' 


Q2n — 


' A 


' 


. A21 


A22 . 





A _ 



\,2mx2n 



where Agt £ C^^" are arbitrary, Q2t is as in Eq.(3.2), and A has the following form 



(3.3) 



^ = ^ ( ylii + A22 ) + ^ ( ^22 - vlii )iei + ^ ( ^21 - ^12 )62 + ^ ( A12 + A21 )ie3 G Q™^'*. (3.4) 
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which can also be stated that for any matrix M G ^2mi<2n^ there must be a unique matrix 
A e Q"*^" such that 

^{A) = M. (3.5) 

Various operation properties on complex representation of complex quaternion matrices can 
easily be derived from Eq.(3.1). 

Theorem 3.2. Let A, B e Q™^", C € Q"^^, and A G C be given. Then 
(a) A = B <^ *(^) = -^{B). 

(h) '^{A + B) = ^{A) + ^{B), ^{AC) = '^{A)^{C), ^ {\A) = ^ {A\) = \^ (A) . 



(c) ^{A*) 



O 



O 



o 



o 



(d) ^'i(A^) = '^*{A), the conjugate transpose of'^{A). 

(e) A = lE2m'^{A)El^, where Eat = [ ( 1 - iei )It, ( 62 + ies )It], t = m, n. 

(f) A is invertible if and only if '^{A) is invertible, in that case, '^[A"^) = '^"^(A) and 

(g) ^{A)ElE2n = El^E2m^{A). 

(h) A is Hermitian if and only if '^{A) is Hermitian over C. 

(i) A is unitary if and only if'^{A) is unitary over C. 

Another universal factorization equality on complex quaternion matrices is established as 
follows. 

Theorem 3.3. Let A = (a^t) G Q"^x", and denote Da = {astl2) e Q2"^x2n y^g^ satisfies 
the following universal factorization equality 



P2mDAP2n 



■0(011) 

V'(ami) 



:= i){A) G C' 



2mx2n 



(3.6) 



where P2t has the form 



P2t = P2t^ = = diag(Q, ■ ■ ■ ,Q), Q = Q-i = Qt = 



1 - iei 62 + ies 
-62 + ie3 1 + iei 



m, n. 



In particular, when m = n, Eq.(3.6) becomes a universal similarity factorization equality over 



Proof. Observe from Eq.(2.1) that Q{astl2)Q = ipidst)- We immediately obtain 

P2mDAP2n = [Q{astl2)Q]mxn = [lp{ast)]mxn = V'(^)- 

which is exactly Eq.(3.6). □ 

The complex matrix tp{A) in Eq.(3.6) is also called the complex representation of A. Clearly 
the two complex matrices "^{A) in Eq.(3.1) and ip(A) Eq.(3.6) are permutationally equivalent, 
that is, there axe two permutation matrices G and H such that G"^{A)H = ip{A). 
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For convenience of application, Eq.(3.6) can be simply stated that for any M G (£^'imx2n^ 
there must exist an A G Q"^x" such that 

V'(^) = M. (3.7) 

Theorem 3.4. Let A, B e Q^^^"^, C G Q^'^'p, XeC be given. Then 

(a) A = B^^{A) =V(5)- 

(b) ^{A + B)=^{A)+^{B), ij{AC) = ij{A)i;{C), ^{XA) = ^{AX) = Xi;{A). 

(c) ll){Ijn)= l2m- 

(d) i^{A^) =ip*{A), the conjugate transpose ofip{A). 

(e) A is invertible if and only if ipiA) is invertible, in that case, ip{A~^) = 'ip~^{A). 

(h) A is Hermitian if and only if ip^A) is Hermitian. 

(i) A is unitary if and only if il}{A) is unitary. 

Just as for complex quaternions, we can define the Moore-Penrose inverse of any m x n 
complex quaternion matrix as follows. 

Definition. Let A G Q^^x" be given. If the following four equations 

AX A = A, XAX = X, (AX)'' = AX, {XA}'' = XA (3.8) 

have a common solution for X, then this solution is called the Moore-Penrose inverse of A, and 
denoted by X = .4+. 

The existence and the uniqueness of the Moore-Penrose inverse of a complex quaternion 
matrix A can be determined by its complex representation ^{A). In fact, according to Theorem 
3.2(a), (c) and (d), the four equations in Eq.(3.8) are equivalent to the following four complex 
matrix equations 

^'(yl)*(X)*(yl) = ^'(^), ^'(X)*(A)^'(X) = *(X), 

[^{A)^{X)]* = ■^{A)^{X), [^{X)^{A)]* = ■^{X)-^{A). 
According to the complex matrix theory, the following four equations 

^'(A)y*(A) = ^'(A), Y^{A)Y = Y, [^(A)F]* = ^(A)F, [Y^{A)]* = Y'^{A) 

have a unique common solution Y = ^^(A), the Moore-Penrose of ^{A). Then by Eq.(3.5), 
there is a unique matrix X over Q such that "^{X) = Y = "^^{A), in which case, this X can be 
expressed as 



X = ^E2nYEl^ = ^E2n^+iA)El^. 



1„ 1 

Correspondingly this X is the unique solution to Eq.(3.8). Hence we have the following. 



Theorem 3.5. Let A G Q^^n given. Then its Moore-Penrose inverse A'^ of A exists 
uniquely, and satisfies the following two equalities 

*(^+) = *+(A), A+ = \E2n^-^iA)El^, 
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where £'2^ = [ ( 1 - iei )It, ( 62 + ie^ )It], t = m, n. 

Through the complex representation in Eq.(3.1), we can define the rank of any A € Q™-^" 
as follows rank(yl) = ^rank[ ^'(A) ]. Obviously the rank of a complex quaternion matrix ^ is a 
fraction if the rank of "^{A) is an odd number. In particular, A G Q"*^"* is invertible if and only 
if rank(74) = m. 

On the basis of the above results we now are able to investigate various kinds of problems 
related to complex quaternion matrices. In the next three sections, we consider three basic 
problems — right eigenvalues and eigenvectors, similarity, as well as determinants of square com- 
plex quaternion matrices. 

4. Right eigenvalues and eigenvectors of complex quaternion ma- 
trices 

As usual, right eigenvalue equation for a complex quaternion matrix A G Q"^" is defined by 

AX = XX, XgQ"^\ AgQ. (4.1) 

If a A G Q and a nonzero X G Q"^-*^ satisfy Eq.(4.1), then A is called a right eigenvalue of A and 
X is called an eigenvector associated with A. In particular, if a A G Q and an X G Q"^^ with 
rank(X) = 1 satisfy Eq.(4.1), then A is called a regular right eigenvalue of A and X is called an 
regular eigenvector associated with A. In this section, we shall prove that any square complex 
quaternion matrix has at least one complex right eigenvalue, and also has at least one regular 
right eigenvalue. To do so, we need some preparation. 

Definition. For any X = Xq + X\ei + ^262 + -^^363 G Q"^^, we call the following complex 
matrix 

:= ^ G C^^^S (4.2) 
uniquely determined by X, the complex adjoint vector of X. 

Lemma 4.1. Let A ^ Q™>^", X G Q"''\ and A G C 6e given. Then a1 = ^{A)t and 

xt = lx. 

Proof. It is easy to see that for any column matrix Y G Q"^^, the corresponding "^{Y) 
and 1^ satisfy the relation = *(y)[ 1, 0]-^. Now applying it to AX and XX we find 

AjI = ^{AX)[1, 0]^ = ^(^)*(X)[1, 0]^ = ^(^)^, 

XX = ^'(XA)[l, 0]^ = ^'(X)A[1, 0]^ = ^A. □ 

Based on the above notation, we can deduce the following several results. 

Theorem 4.2. Let A G Q"^" be given. Then all the eigenvalues of'^{A) are right eigenvalues 
of A, and all the eigenvectors of'^{A) can be used for constructing eigenvectors of A. 

Proof. Assume that A G C and O 7^ Y G C^"^^ satisfy 

*(A)y = YX. 



Xq + X^i 
X2 - Xsi 



9 



Then we set X = -E'2n^ £ Q"^-^, where i?2n = [(1 ~ 'iei (e2 + ies)In], and it satisfies 
E2nE2n = Aln- Combining the above results with Theorem 3.2(e) and (g), we obtain 

AX = ^E2n^{A)El^E2nY = ^E2n4^E2n^{A)Y = E2nY X = XX, 

which shows that A is a right eigenvalue of A and X = E2nY is an eigenvector of A associated 
with this A. So the conclusion of the theorem is true. □ 

Conversely, we have the following result. 

Theorem 4.3. Let A G Q"^" be given. Then all the complex right eigenvalues of A are 
eigenvalues of'^{A), too. 

Proof. Assume that A € C and O ^ X £ Q"^^ satisfy AX = XX. Then applying Lemma 4.1 
to the both sides of this equality, we find '^{A)jt = jtx, which shows that A and jt are a pair 
of eigenvalue and eigenvector of '^{A). □ 

Next are two results on the regular right eigenvalues and eigenvectors of complex quaternion 
matrices. 

Theorem 4.4. Let A G Q"^" be given. Then A has at least one regular right eigenvalue. 

Proof. We first assume that the complex representation ^(A) has at least two linearly 
independent eigenvectors Yi and I2, i-©-, 

*(A)yi = yiAi, ^iA)Y2 = Y2X2, (4.3) 

where rank[yi, Y2] = 2, Ai and A2 are two complex numbers with Ai = A2, or Ai 7^ A2. Then 



^iA)[Yi, Y2] = [Fi, Y2 



Ai 
A2 



According to Eq.(3.5), there must exist an X G Q"^-*^ such that 

^'(X) = [Yi, Y2], rank(X) = ^rank*(X) = 1, 



meanwhile there must exist a A G Q such that 

^(A) = 



Ai 
A2 



(4.4) 



(4.5) 



(4.6) 



Putting Eqs.(4.5) and (4.6) into Eq.(4.4) and then applying Theorem 3.2(a) to it, we obtain 

^(A)^(X) = ^(X)^'(A) =^AX = XX, and rank(X) = 1, 

which shows that A and X are a pair of regular right eigenvalue and eigenvector of A. 

Next assume that the complex matrix ^(A) has only one eigenvalue Ai and only one linearly 
independent eigenvector Yi associated with Ai. Then according to complex matrix theory, there 
exists another vector Y2 over C such that 



^{A)[Yi, Y2] = [Yi, Y2] 



Ai 1 
Ai 



(4.7) 
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where rank[Yi, I2] = 2. Then it follows by Eq.(3.5) that there must exist a A G Q and an 
X G Q"^^ such that 



*(A) 



Ai 1 
Ai 



, ^{X) = [Yi, Y2], rank(X) = -rank*(X) = 1. 



(4.8) 



In that case, combining Eqs.(4.7) and (4.8) and Theorem 3.2(a), we get 

^{A)'^{X) = ^'(X)*(A) ^ AX = XX, r{X) = 1, 
which shows that A and X are a pair of regular right eigenvalue and eigenvector of A. 



□ 



Theorem 4.5. Let A G Q"^" be given. Then from any regular right eigenvalue of A we can 
derive eigenvalues of'^{A). 

Proof. Assume that A = Aq + Aici + A2e2 + Xse^ e Q and X e Q"^^ with rank(X) = 1 satisfy 

AX = XX. (4.9) 
If A G C in Eq.(4.9), then from Theorem 3.3(a) and (b) we get 



*(y4)*(X) = *(X)*(A) = -^{X) 



X 
A 



which clearly shows that A is an eigenvalue of ^{A), and ^{X) arc two eigenvectors of ^(A). If 
A ^ C in Eq.(4.9), then by Theorem 2.6 we know that this A can be expressed as 



A = p[Ao + r(A)ei ]p-\ 

when T^{X) = Af + A| + A| 7^ 0, or 

A = g ^ Ao - + ^ies ^ q~^, 



(4.10) 



(4.11) 



when Af + A| + A| = 0. Under the condition in (4.10), the equality in Eq.(4.9) can be equivalently 
expressed as 

A{Xp) = (Xp)[Ao + r(A)ei]. (4.12) 
Applying Theorem 3.2(b) to the both sides of Eq.(4.12) yields 



^{A)<if{Xp) = ^i{Xp)^Xo + r(A)ei ] = ^{Xp) 



Ao + r(A)i 

Ao - T{X)i 



which shows that Aq ± ^(A)^ are two eigenvalues of '^{A), and ^{Xp) are two eigenvectors of 
'^{A). Under the condition in Eq.(4.11), the equality in (4.9) can equivalently be expressed as 



A{Xq) = {Xq)(^Xo-^e2 + ^iesy 
Then applying Theorem 3.2(b) to the both sides of Eq.(4.13) yields 
^{A)^iXq) = ^{Xq)^ (Xo- \ei + ^^63] = *(Xg) 



(4.13) 



Ao 1 
Ao 
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which shows that Aq is an eigenvalue of "^{A), and the first column of "^{Xp) is an eigenvector 
of*(A). □ 

5. Similarity of complex quaternion matrices 

Two square matrices A and B of the same size over Q arc said to be similar if there exists 
an invertible matrix X over Q such that X~^AX = B. Based on the results in the preceding 
two subsections, we can easily find a simple result charactering the similarity of two complex 
quaternion matrices. 

Theorem 5.1. Let A, B & Q"^'^ be given. Then the following three statements are equivalent: 

(a) A and B are similar over Q. 

(b) ^{A) and ^{B) are similar over C 

(c) ipiA) and ip{B) are similar over C. 

Proof. Assume first that A ^ B over Q. Then there is an invertible matrix X such that 
X-^AX = B. Now applying Theorem 3.2(b) and (f) to its both sides yields ■^~^{X)^{A)^(X) = 
^{B), which shows that ~ over C. Conversely, assume that ^{A) ~ *(-B) over 

C. Then there is an invertible matrix Y G C^'''ix2n ^^^^^i that Y-'^'^{A)Y = For this 

Y, by (3.5) there must be an invertible matrix X G Q"^" such that "^{X) = Y. Thus, 
^(X-^)^(A)^{X) = ^'(S), and consequently X'^AX = B, which shows that A B. The 
equivalence of (a) and (c) can also be shown in the same manner. □ 

Based on this result, we can extend various results on similarity of complex matrices to 
complex quaternion matrices. 

Theorem 5.2. Let A G Q"^"^ be given. Then A is similar to a diagonal matrix over Q if and 
only if the sizes of the Jordan blocks in the Jordan canonical form of the complex representation 
ip{A) of A are not greater than 2. 

Proof. Assume that A is diagonalizable over Q, i.e., there is an invertible matrix P over Q 
such that 

P-^^P = diag(Ai, A2, A„), (5.1) 
where Ai — A„ G Q. Then by Theorem 5.1, we obtain 

V'-^(P)V'(^)V'(^) =diag(V(Ai), V'(A2), V'(An)), (5.2) 

where V'(At) is a 2 x 2 complex matrix. This equality implies that the sizes of Jordan block in 
the Jordan canonical form of V'(^) are not greater than 2. Conversely assume that there is an 
invertible matrix G over C such that 

G~V(^)G = diag( J(|Ui), J^Hk), fJ-k+i, Mr), (5.3) 

where J{iit){t = 1 — k) is a 2 x 2 complex Jordan block, jik+i — A*r are complex eigenvalues of 
V'(A) with 2k + {r — k) = 2n. In that case, by Eq.(3.7) we know that there is an invertible 
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matrix P over Q such that (j){P) = G, there is a Xt over Q such that ip{\t) = J{f^t){t = 1 — fc)) 
and there are A^+i, • • • , A,- € Q such that 







^(A 



A:+2j 



^^k+A 



fJ.r-1 
Hr 



(5.4) 



In that case, let D = diag( Ai, A2, • • • , A„ ). Then Eq.(5.3) becomes 



Hence A is diagonahzable over Q. □. 



P-^AP = D. 



The result in the above theorem alternatively implies that if the size of a Jordan block in 
the Jordan canonical form of the complex representation ip(^A) of A is greater than 2, then A 
can not be diagonahzable over Q. 

Theorem 5.3. Let A E Q*^^" be given. Then A is similar to a complex matrix over Q if and 
only if the complex representation ^{A) of A is similar to a block diagonal matrix diag( J, J) 
over C, where J £ C""^". 

Proof. If yl ~ J G C"''", then it follows by Theorem 5.1 and Eq.(3.1) that ^{A) ~ ^{J) = 
diag( J, J). Conversely if ^'(A) ~ diag( J, J) with J E C"^", then it follows by (3.1) that 
^'(J) = diag( J, J). Consequently A^ J follows by Theorem 5.1. □ 



6. Determinants of complex quaternion matrices 



As one the most fundamental problems in the theory of quaternion matrices, the determi- 
nants of complex quaternion matrices, including the determinants of real quaternion matrices, 
have been considered by lots of authors from different aspects, see, e.g., [^, |^], P and [p^ . 
Among all of them, a direct and simple method for defining determinants of quaternion matrices 
is through their representations in the central fields of the corresponding quaternion algebras. 
In [|^, Zhang presents this kind of definition for determinants of real quaternion matrices, and 
demonstrates the consistency of his definition with some other classic definitions on determinants 
based on the products and sums of entries in matrices. 

Following the introduction of the universal similarity equalities for generalized quaternion 
matrices, we shall easily find that it is a most reasonable method to define determinants of 
quaternion matrices through their representations in the central field of corresponding quaternion 
algebra. In fact, from Eq.(3.1), we know that for any A € Q"^" 

Q2ndiag( A, A )Q^^ = ^{A) G C^-x^". (6.1) 

Therefore if we hope that determinants of quaternion matrices over Q satisfy the following two 
basic properties 

det(MiV) = det MdetA^ and det/„ = 1 

for any M, N € Q"^", then as a natural consequence of these two properties, the determinant 
of the diagonal matrix diag(y4. A) satisfies the following equality 

det[diag( A A)]= det[Q2ndiag(A, A)Q:^^] = det^'(A). (6.2) 
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From this equality we see that there seem only two kinds of natural choices for the definition of 
determinant oi A E Q"^", namely, define 



detA := l^-l^)!, 



(6.3) 



or 



detA:=\'^{A)\2, (6.4) 

where |^'(A)| is the ordinary determinant of the complex matrix ^{A). For simplicity, we prefer 
Eq.(6.3) as the definition of determinant of matrix over Q and call it the central determinant of 
matrix A, and denote it by \A\c. 

Some basic operation properties on the central determinants of complex quaternion matrices 
are listed below without proofs. 

Theorem 6.1. Let A, B £ Q"^", X e C, and fj, e Q be given. 

(a) //A € C"^", i/ien = l^p. 

(b) A is invertible <^=^ \A\c 7^ 0. 

(c) \ABU=\AU\BU. 



(d) \XA\, = \AX\, -- 

(e) IMIc = = 

(f) \A-%=JA\-\ 

(g) l^^lc = |^le- 
an 

(h) If A 



X'^"\A\c. 

n?"'{lJ,)\A\c, where n{ii) is the weak norm of ji. 



(i) If A 



Ai * 
A2 

(j) IfA^B, then \A\c = \B 



, then \ A\c = n(aii)n(a22) • ■ ■ n(a„„). 



, where Ai and A2 are square, then \A\c = \A1\c\A2\ 



By means of the determinants of complex quaternion matrices defined above, we can easily 
define the central characteristic polynomial of any A G Q'*^" as follows 

PA{X) = \Xl2n-'^{A)\, (6.5) 

which is a complex polynomial of degree 2n. From it we easily get the following. 

Theorem 6.2(Cayley-Hamilton theorem). Let A G Q"^" be given. Then pa{A) = 0. 

Proof. Since pa(A) = \XI2n - *(^)|, so pa[^{A)] = 0. Putting Eq.(6.1) in it and simplifying 
the equality yields 

PA[diag{A, A)] = diag(pyi(^), Pa{A) ) = 0, 
which implies that pa{A) =0. □ 



7. Conclusions 
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In the article, we have estabhshed a fundamental universal similarity factorization equality 
over the complex quaternion algebra Q. This equality clearly reveals the intrinsic relationship 
between the complex quaternion algebra Q and the 2x2 total matrix algebra, and could serve 
as a powerful tool for investigating various problems related to complex quaternions and their 
applications. In addition to the results in Sections 3 — 6, we can also apply Eqs.(2.4), (3.1) 
and (3.6) to investigate some other basic topics related to complex quaternion matrices, such 
as, singular value decompositions, norms, numerical ranges of complex quaternion matrices and 
so on. Finally we point out that the equality (2.4) can also extended to the complex Clifford 
algebra C^, and a set of perfect matrix representation theory on the complex Clifford algebra 
Cn can explicitly established. We shall examine this topic in another paper. 
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